Critical Temperature of Weakly Interacting Dipolar Condensates 
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We calculate perturbatively the effect of a dipolar interaction upon the Bose- Einstein condensation 
temperature. This dipolar shift depends on the angle between the symmetry axes of the trap and the 
aligned atomic dipole moments, and is extremal for parallel or orthogonal orientations, respectively. 
The difference of both critical temperatures exhibits most clearly the dipole-dipole interaction and 
can be enhanced by increasing both the number of atoms and the anisotropy of the trap. Applying 
our results to chromium atoms, which have a large magnetic dipole moment, shows that this dipolar 
shift of the critical temperature could be measured in the ongoing Stuttgart experiment. 

PACS numbers: 03.75.Hh, 31.15.Gy, 51.30.+i 



Ultracold atomic quantum gases are many-body sys- 
tems in which macroscopic quantum phenomena can be 
studied experimentally over a wide range of controllable 
interactions [1-3]. For the original alkali atomic Bosc- 
Einstein condensates (BECs) it has been sufficient to 
describe the dominant two-particle interaction by a lo- 
cal isotropic contact potential. Recently, a new type of 
nonlocal anisotropic interaction has been made accessible 
to detailed study by the formation of a BEC in a dipo- 
lar quantum gas of 52 Cr [4]. Chromium atoms possess 
a magnetic dipole moment of 6 to_b, where is the 
Bohr magneton, i.e. it is around six times larger than 
those of alkali atoms. It stands for a whole class of high 
spin atoms like rare earth atoms with large magnetic 
dipole moments (Dy [10 TOb], Ho [9 tub], Eu [7 tub], Tb 
[10 m B ], Er [7 m B ], Mo [6 m B ], Mn [5 m B ], Tm [4 m B ], 
Pr [3.3 tub}). Note that these elements have been already 
cooled by buffer gas and evaporative cooling techniques 
down to mK temperatures [5]. Besides that Er [7 m B ] 
has recently been laser cooled [6] . Other many-body sys- 
tems with dipolar interactions are, for instance, Rydberg 
atoms [7, 8] or atomic condensates where a strong electric 
field induces electric dipole moments of the order of 10 -2 
Debye [9]. Body-centered dipole moments in heteronu- 
clear molecules are much larger with typical values of 1 
Debye, so their dipolar effects could be a few hundred 
times stronger than those of chromium atoms [10, 11]. 
Such a gas of ultracold heteronuclear molecules is pro- 
duced either by sophisticated cooling and trapping tech- 
niques [12, 13] or by photoassociation [14-16]. For all 
those dipolar systems the total two-particle interaction 
potential is modelled by 
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mass. A more general pseudopotential for anisotropic in- 
teractions has recently been introduced in Refs. [17, 18] 
which is nonlocal in momentum space. Since there ex- 
ists up to now no experimental evidence for any dipolar 
shape resonance, where such a more general pseudopo- 
tential could be relevant, our model (1) is valid for all 
current and many future experimental situations. 
The dimensionless measure of the strength of the dipole- 
dipole interaction with respect to the s-wave scattering 
is e DD = iiom 2 M/(12iTh 2 a) [19]. For the 52 Cr conden- 
sate it has the value €dd — 0.144 [20], so the magnetic 
dipolar interaction represents only a small correction to 
the contact interaction. Nevertheless, it has interesting 
consequences due to its anisotropy, as has been observed 
in a recent expansion experiment [21]. Furthermore, the 
magnetic dipole-dipole interaction can be varied within a 
limited range with the help of rotating magnetic fields as 
proposed by Ref . [22] . Combining this technique with the 
now known 14 Fcshbach resonances of chromium atoms 
[20] will allow experiments where the interaction varies 
from only contact to purely dipolar. In this way, many 
interesting predicted dipolar phenomena should be ob- 
served. Among them are, for instance, the stability of 
the ground state of a dipolar BEC [23-25] or its excita- 
tion spectrum [26]. Note that it has been recently sug- 
gested in Refs. [27, 28] that the s-wave scattering length 
a could strongly depend on the dipole moment. However, 
there it has also been shown that for dipolar interaction 
strengths, which are not larger than the s-wave scattering 
strength, the latter is only rescaled by a moderate factor 
and remains positive. For the calculations in the present 
work we assume that this condition is fulfilled. 
In this letter, we investigate how the critical tempera- 
ture of a dipolar BEC depends on the dipole-dipole in- 
teraction. Consider an atomic gas trapped in a cylinder- 
symmetric harmonic potential 



Here m denotes the magnetic (electric) dipole moment of 
the atoms, /Uo stands for the magnetic field constant (the 
reciprocal electric field constant), and M is the atomic 
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whose dipole moments m have the angle a with the pr- 
axis, i.e. m = to (sin a, 0, cos a). In particular, we are 
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FIG. 1: Symmetry axes of harmonic trap (2) and two-particle 
interaction (1) are parallel and perpendicular in configuration 
I and II. 



interested in the two extreme configurations I and II 
where the symmetry axis of the dipole moments is paral- 
lel (a = 0) and perpendicular (a = ir/2) to the symmetry 
axis of the harmonic trap (2), respectively. Figure 1 il- 
lustrates the two cases for uu < uij_, where the atomic 
dipole moments lead to a residual attractive and repul- 
sive interaction, respectively. At first we show that the 
corresponding critical temperature in configuration I and 
II is shifted above and below the value of a pure contact 
interaction. Subsequently we suggest to determine the 
difference of the critical temperatures in both configura- 
tions I and II to cancel out the influence of the isotropic 
contact interaction. We shall estimate which experimen- 
tal parameters allow us to enhance this signal. For a pure 
contact interaction, the harmonic trap suppresses long- 
wavelength fluctuations, so the leading shift of the crit- 
ical temperature can be calculated perturbatively (see, 
for instance, the recent work [29] and the references cited 
therein) . We expect that this reasoning also holds for our 
model interaction (1). Although the dipolar interaction 
is nonlocal, its scaling properties are the same as for a 
contact interaction. 

In this letter we apply Feynman's diagrammatic tech- 
nique of many-body theory [30, 31] and expand the 
grand-canonical free energy as 
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The first term denotes the contribution of an ideal Bose 
gas at temperature T — 1/keP where the harmonic trap 
potential (2) is treated semiclassically [32] 
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It contains the chemical potential /i, the geometric mean 
of the trap frequencies ui — (ojii WjJ 1 / 3 , and the polyloga- 
rithmic function C, a (z) — Ysn=i z " / n °- The two diagrams 
in Eq. (3) represent the direct and the exchange vacuum 
contribution, respectively, and have to be evaluated ac- 
cording to the Feynman rules: a straight line with an 
arrow represents the semiclassical interaction-free corre- 
lation function: 
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where io m = 2-Km/Ti{3 denotes the Matsubara frequency. 
The two-particle interaction potential enters via the dia- 
gram 
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The grand-canonical free energy (3) is studied as a func- 
tion of temperature T for fixed particle number N = 

The phase transition, where a macroscopic occupation 
of the ground state sets in, occurs when the correlation 
function of the system diverges. From Eq. (5) we can see 

that this happens in the interaction- free case at /j,^ = 0, 
as the divergence appears at the minimum of the trap po- 
tential V(x) for vanishing Matsubara frequency tu m and 
momentum p. In the presence of a 2-particle interaction, 
the full correlation function follows from a formula sim- 
ilar to (5) where the chemical potential \i is shifted by 
the self-energy: p — ► /i + 7i£(p, w m ;x). It is defined by 
the Fourier-Matsubara transformed 
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The critical chemical potential reads now [33] 

V(x) - ftE(p,0;x) 
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which leads to /j c = — ft, 11(0, 0; 0) up to first order. 
Thus, evaluating the particle number N — N(ji) at 
the critical chemical potential fi = fi c yields the follow- 
ing first-order shift of the critical temperature with re- 
spect to the interaction-free critical temperature Tj — 
TiCo[N/aZ)] 1/3 /kB: 

[3cos a - 1] / ( -JL ) - .(10) 
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Here £(a) = Yl^L^ l/n a is the Riemann zeta function 
and A = (2-rrft 2 0/M) 1 ' 2 the thermodynamic de Broglie 
wave length. The dimensionless prefactor c$ for the 5- 
interaction has the value 
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with the generalized Riemann zeta function £(a,b,c) — 

J2™=i En/=i l/n a n' ( n + n ') C i whereas the dimension- 
less prefactor / (k) for the dipole-dipole interaction is a 
function of the ratio of the trap frequencies k = lu»/uj±: 
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FIG. 2: a) Shift of the critical temperature AT C with respect 
to the interaction- free critical temperature ri ' for a 52 Cr gas 
in a harmonic trap with frequencies ui^ — 27T-138 Hz and ui± = 
2-7T • 486 Hz: without (straight line) and with (dashed lines) 
magnetic dipole-dipole interaction for the configurations I and 
II of Figure 1. b) Difference of the temperature shifts increases 
with particle number N and geometric mean frequency Q. 
The respective dots indicate the parameters of the Stuttgart 
52 Cr experiment. 
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FIG. 3: a) Shift of the critical temperature AT C with respect 
to the interaction-free critical temperature Tj ' for N = 10 5 
52 Cr-atoms gas versus anisotropy parameter wi\/lj± of the 
harmonic trap without (straight line) and with (dashed lines) 
magnetic dipole-dipole interaction for the configurations I and 
II of Figure 1. b) Difference of the temperature shifts versus 
anisotropy parameter uj\\/u±_ for CjN 1 '' 3 = 0.93 • 10 s Hz (solid 
line) and OjN 1 ^ 3 = 3-10 5 Hz (dashed line). The respective dots 
indicate the parameters of the Stuttgart 52 Cr experiment. 



This function was already used in Ref. [19] to describe 
the mean-field magnetic dipole-dipole energy for a cylin- 
drically symmetric BEC. Note that /(k) tends asymp- 
totically to — 2 for k — ► oo and 1 for k — ► 0, respectively. 
The general physical implications of our first-order per- 
turbative result (10) are as follows. Without the dipole- 
dipole interaction, i.e. at m = 0, the critical tempera- 
ture is shifted downwards with the dimensionless prefac- 
tor (11). This result for the isotropic contact interaction 
was originally derived in Ref. [34] within a mean-field 
approach and confirmed experimentally by investigating 
the onset of Bose-Einstein condensation in the hypcrfinc 
ground state of 87 Rb [35]. Our result (10) for m ^ 
shows that the critical temperature is increased in con- 
figuration I (a = 0) twice as much as decreased in config- 
uration II (a = 7r/2) of Figure 1 for the trap anisotropy 
wii < uix, due to the dipole-dipole interaction. Of course, 
changing the angle a allows us to tune the dipolar effect 
between these maximal and minimal values. In partic- 
ular, we read off from Eq. (10) that the dipolar shift 
of the critical temperature vanishes for the magic angle 
a = arccos(l/\/3) = 54.7° [22]. 

Now we discuss the consequences of our results for the 
ongoing experiments on the Bose-Einstein condensation 
of 52 Cr-atoms at the University of Stuttgart, where the 
trap frequencies are u)\ — 2tt ■ 581 Hz, 102 — 2tt • 406 Hz, 
CJ3 = 2ir ■ 138 Hz, so that the geometric mean frequency 
is u> — 2tt ■ 319 Hz. The total number of atoms is 
N = 100 000, yielding an interaction-free critical tem- 



perature of about T, 



(0) 



670 nK. The corresponding 



finite-size correction was calculated in Refs. [36, 37] 
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where to is the arithmetic mean of the trap frequencies 
ZJ w (u>i + ll>2 + /3 = 2tt ■ 375 Hz. Thus, the finite-size 
correction of Tc°^ in the Stuttgart experiment amounts 
to —1.8 %. This is to be compared with a shift of the 
critical temperature due to the contact and the magnetic 
dipole-dipole interaction following from formula (10). 



The s-wave scattering length of the 52 Cr-atoms is 
a = 105 as [20], i.e. roughly 2 % of the thermodynamic 

de Broglie wave length = 5 598 as- The correspond- 
ing downwards shift of the critical temperature amounts 
then to 6.4 %. This is modified by the dipole-dipole 
interaction depending on the experimental set-up. In 
the Stuttgart experiment, the trap potential has actu- 
ally three different frequencies but, since u>\ and 0J2 are 
not far apart, we may identify lo±_ = ^JuJ\uj2 = 2ix ■ 486 Hz 
and wii = u>3, yielding ll>\\/uj± — 0.284. In configuration I 
we obtain from (12) the result f(uin/ui±) = 0.733, which 
leads via Eq. (10) to an extra upward shift of the critical 
temperature by 0.34 % due to the magnetic dipole-dipole 
interaction. The downward shift in configuration II is 
half as big. Figure 2a plots the resulting total shift of the 
critical temperature AT C for the 52 Cr gas with respect to 
the interaction-free critical temperature versus the 
particle number N. Both the finite-size corrections and 
the contact interaction lead to the main shift, on top of 
which the small dipolar effect is seen. Figure 3a shows 
how the same shifts depend for N = 10 chromium atoms 
on the anisotropy parameter u)\\/w_\_. The directions of 
the shifts change sign at the isotropy point u)\\ = wi. 
The above results suggest to plot the difference between 
the critical temperatures of the two configurations I and 
II. This eliminates the isotropic effects of both the con- 
tact interaction and the finite-size correction, thus it ex- 
hibits most clearly the magnetic dipole-dipole interac- 
tion. For N = 10 5 and Co = 2n ■ 319 Hz, the difference 
amounts to a net effect of 0.51 % of the interaction-free 
critical temperature T c . One possibility to enhance the 
difference of the critical temperatures is to choose a con- 
venient anisotropy strength u>»/u>x of the harmonic trap 
potential as seen in Figure 3b. Furthermore, this dipo- 
lar effect increases with the total atom number TV and 
the geometric mean frequency u) as shown in Figure 2b, 
which therefore needs to be as large as possible. An ex- 
perimentally feasible increase of the particle density in 
the trap could lead to a dipolar effect of more than 2 %, 
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see dashed curve of Figure 3b. 

At present, the best experiments which measure the crit- 
ical temperature of a Bose-Einstein condensate involve 
error bars of 5 % which represent the total systematic 
and statistical errors [4, 35]. The systematic errors, how- 
ever, can be eliminated by our suggestion to measure the 
difference of two critical temperatures. We expect that 
the remaining statistical errors can be reduced to 1 % by 
averaging both atom number and critical temperature 
over many measurements. 

So far, the dipolar nature of the chromium BEC has 



only been resolved in expansion experiments [21]. The 
analysis of this letter shows that it could be possible to 
detect a signal of the underlying magnetic dipole-dipolc 
interaction also by measuring the critical temperature. 
Furthermore, our results will be useful for other dipolar 
systems with a tunable dipole moment, like heteronu- 
clear molecules in low vibrational states [12-16], where 
the dipolar effect will be larger. 

We thank A. Berra and S. Kling for discussions and DFG 
Priority Program SPP 1116 as well as SFB/TR 21. 
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